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Abstract. We show that a first order perturbation A{x) ■ D + q{x) of the 
polyharmonic operator (—A)™, m > 2, can be determined uniquely from the 
set of the Cauchy data for the perturbed polyharmonic operator on a bounded 
domain in M", n > 3. Notice that the corresponding result does not hold in 
general when m = 1. 



1. Introduction and statement of results 

Let r2 C M", n > 3, be a bounded domain with C°°-boundary, and consider the 
polyharmonic operator (—A)*", m > 1, which is a positive self-adjoint operator 
on L'^iyi) with the domain 

n ifo"(^)> Ki^) = {ue H"'{n) : = 0}. 

Here 

7M = {u\dn, di,u\on, d'^^'^ulan) 

is the Dirichlet trace of u, u is the unit outer normal to the boundary dfl, and 
H^{Q) is the standard Sobolev space on Q, s E M. 

Consider the operator 

n 

i=i 

with A = iAj)i<j<n e W^'°°{nX") and q G L°°{n,C). Here D = r^V. Viewed 
as an unbounded operator on L^(f2) and equipped with the domain H'^"^{Q) fl 
H^{fl), the operator is closed with purely discrete spectrum. 

Let us make the following assumption, 

(A) is not an eigenvalue of CaA^, D) : H'^^'^n) n Hl^{n) ^^(fi). 

Under the assumption (A), for any / = (/o, /i, • ■ ■ , /m-i) ^ 'HP''^^^{dVi) : = 
nr=V ^^'""^"^/^(<9fi), the Dirichlet problem 

Ca = in VL, 

(1.1) 

'ju = f on dQ, 
1 
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has a unique solution u G H'^"^{Q). Introducing the Neumann trace operator 7 
by 

2m-l 
j=m 

7«=(5>lac,...,5^"^«|an), 
see [IB] , we define the Dirichlet-to-Neumann map MA,q = ■^A,q by 

A^,, : •H°''"-^(9r]) ^ Ar4,,(/) = 7M, 

where u G H'^"^{Q) is the solution to the Dirichlet problem (11.11) . For future 
references, let us also introduce the set of the Cauchy data CA,q = C^^^ for the 
operator CA,q defined as follows, 

C^_g = {(7m,7m) : u G H^'^in), CA,qU = in n}. 

When the assumption (A) holds, the set CA,q is the graph of the Dirichlet-to- 
Neumann map MA,q- 

In this paper we are concerned with the inverse problem of recovering the first 
order perturbation A{x) ■ D + q{x) in Q from the knowledge of the Dirichlet-to- 
Neumann map MA,q on the boundary of Q. 

As it was noticed in [38], when m = 1, we have 

e-''l'CA,qe'^ = C^^g, e'^^A,q(^-'^ = ^A,q - td.^P, 

A = A + 2V'iIj, q = q + A-V^+\Vip\'^ -iA-^. (1.2) 

Thus, N'A^q = -^Aq "when ip G C^(f2) is such that ilj\dn = <9,^'0|9n = 0. We may 
therefore hope to recover the coefficients A and q from boundary measurements 
only modulo a gauge transformation in (II ■2p . 

Starting with the pioneering paper [38j, inverse boundary value problems for 
first order perturbations of the Laplacian have been extensively studied, usually 
in the context of magnetic Schrodinger operators. In [38] it was shown that 
the hope mentioned above is justified, and the magnetic field and the electric 
potential are uniquely determined by the Dirichlet-to-Neumann map, provided 
that the magnetic field is small in a suitable sense. The smallness condition was 
removed in [33] in the case of C°°-smooth magnetic and electric potentials, and 
also for C^-compactly supported magnetic and L°° electric potentials, see also 
[B] . The regularity assumptions on the potentials were subsequently weakened in 
the works [371 SI] ■ The case of partial boundary measurements for the magnetic 
Schrodinger operator was studied in [TJ [25] . 

The purpose of this paper is to show that the obstruction to uniqueness coming 
from the gauge invariance (II. 2p when m = 1, can be eliminated by passing to 
operators of higher order, and the unique determination of the coefficients A and 
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q becomes then possible. Throughout the paper we shall assume therefore that 
m > 2. Our first result is as follows. 

Theorem 1.1. Let Q C M", n > 3, be a bounded domain with C°° -boundary, 
and let AW,^^^) g H^1'°°(M", C") n ^(n, C") and gW,g(2) ^ L~(fi,C), be such 
that the assumption (A) is satisfied for both operators. If Afj^w = A/'^(2) ,^(2), 
then A(^) = ^(2) and q^^^ = q^^^ in Vt. 

The proof of Theorem 11.11 exploits complex geometric optics solutions for the 
equations = 0, j = 1,2. Following [71 [21], the construction of such 

solutions will be carried out using Carleman estimates. We remark here that 
starting with the fundamental paper |40j, complex geometric optics solutions have 
been a central tool in establishing uniqueness results in elliptic inverse boundary- 
value problems. 

Dropping the assumption that A^^'^ = 0, j = 1,2, along dVL in Theorem II. H we 
have the following result. 

Theorem 1.2. Let C M", n > 3, &e a bounded domain with C°° -boundary, and 
let G C°°(n,C") and q^^\q^'^^ G C°°(n,C), be such that the assumption 

(A) is satisfied for both operators. If MAW^qW = J^Ai^) ,q(^') > then A^^^ = A^'^'> and 

The idea of the proof of Theorem 11.21 is to reduce it to Theorem 11.11 by using the 
boundary reconstruction of the vector field part of the perturbation. When doing 
the boundary reconstruction, similarly to [231 ISll |33] , we use pseudodifferential 
techniques, which motivates the need to require C°° smoothness assumptions on 
the coefficients of the perturbed operator. 

Finally, it may be interesting to notice that the boundary reconstruction becomes 
unnecessary when the boundary of the domain Q is connected. We have the 
following result. 

Theorem 1.3. Let Q C M", n > 3, be a bounded domain with connected C°°- 
boundary, and letA^^\A^^^ G iyi'~(fi,C") andq^^\q^^^ G L°°(fi,C), be such that 
the assumption (A) is satisfied for both operators. If N^w^qW = A/'^(2) g(2), then 

^(1) _ ^(2) ^(1) _ ^(2) ^_ 

The inverse boundary value problem of the recovery of a zeroth order perturbation 
of the biharmonic operator has been studied in [201 122] , and unique identifiability 
results were obtained, similarly to the case of the Schrodinger operator. The 
areas of physics and geometry where higher order operators occur, include the 
study of the Kirchhoff plate equation in the theory of elasticity, and the study of 
the Paneitz-Branson operator in conformal geometry, see [10]. Inverse boundary 
value problems for differential perturbations of the biharmonic, or more generally. 
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polyharmonic, operator are therefore natural to consider, and the present paper 
is meant as a step in this direction. 

We would also like to mention the important results of [8], concerned with the 
inverse boundary value problem for a first order matrix perturbation of the Lapla- 
cian, considered in a smoothly bounded convex domain in M", n > 3. The results 
of [8] show that it is possible to recover the matrix potentials up to a gauge trans- 
formation, given by a smooth invertible matrix. In our study of the perturbed 
polyharmonic operator, rather than reducing it to a system, we adopt a direct 
approach, which has the merit that the issue of the gauge equivalence does not 
arise. 

Finally, we would like to point out that the results obtained in this paper can 
be viewed as generalizations of the corresponding results for second order equa- 
tions, encountered in electrical impedance tomography, see [21 El |32l [39] for the 
two dimensional case, and IHl |311 [Ml [10] for the higher dimensions, as well 
as in inverse boundary value problems and inverse scattering problems for the 
Schrodinger equation |H [Til HU [321 [311 HO] , and in elliptic inverse problems on 
Riemannian manifolds, [T71 [TH |2Z1 [281 129]. There are also counterexamples for 
uniqueness of inverse problems with very non-regular electric fields [12] and mag- 
netic fields [15]. These counterexamples are closely related to the so-called invis- 
ibility cloaking, see [III [131 [H |26l [SOl [35] . 

The plan of the paper is as follows. In Section 2 we construct complex geometric 
optics solutions of the equation CA,qU = in fi, which are instrumental in proving 
Theorems I1.1H1.3I The proof of Theorem 11.11 is then given in Section 3, while 
the boundary reconstruction of the vector field part of the perturbation and the 
proof of Theorem 11.21 are the subjects of Section 4. The proof of Theorem 11.31 
is finally given in Section 5. Appendix A contains a characterization of curl-free 
vector fields which may be of some independent interest. 

2. Construction of complex geometric optics solutions 

Let VL C M*^, n > 3, be a bounded domain with C°°-boundary. Following [71 [21], 
we shall use the method of Carleman estimates to construct complex geometric 
optics solutions for the equation CA,qU = in with A G VF^'°°(fi,C") and 
gGL°°(n,C). 

First we shall derive a Carleman estimate for the semiclassical polyharmonic op- 
erator (— /i^A)™, where > is a small parameter, by iterating the corresponding 
Carleman estimate for the semiclassical Laplacian — /i^A, which we now proceed 
to recall following [21]. Let VL be an open set in such that Vt dd Vt and 
e C°°(n,M). Consider the conjugated operator 
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and its semiclassical principal symbol 

Following [24], we say that 93 is a limiting Carleman weight for —h?/S. in f2, if 
V(/5 7^ in f2 and the Poisson bracket of Rep,^ and Imp^ satisfies, 

{Repy,,Imp^}(a;,^) = when j9^(x,^)=0, (x,^)Gfix]R". 

Examples are linear weights ^{x) = a-x, a G M", |a| = 1, and logarithmic weights 
<f{x) = log |x — xol, with xo ^ fi. In this paper we shall only be concerned with 
the case of linear weights. 

In what follows we shall equip the standard Sobolev space i/'^(M"), s G M, with 
the semiclassical norm HmUh^^j = \\{hD)'^u\\i2. Here (^) = (1 + 1^^)"'^''^. We shall 
need the following result, obtained in [21]. 

Proposition 2.1. Let be a limiting Carleman weight for the semiclassical 
Laplacian on Q. Then the Carleman estimate 

h 

\\e^{-h^A)e'^u\\H^ > ^\\u\\j,s+., (2.1) 

Cs,n ^'^i 

holds for all u G C^{VL), s G M and all h > small enough. 

Iterating the Carleman estimate (12.11) m times, m > 2, we get the following 
Carleman estimate for the polyharmonic operator, 

urn 

\\e^{-h^Are-^u\\H^ > ■^\\u\\rrs+^, (2.2) 

for all u G C^{Q), s G M and h > small. Since we are dealing with first order 
perturbations of the polyharmonic operator, the following weakened version of 
(12. 2p will be sufficient for our purposes, 

um 

\\eH-h'Ay-e-U\\H:^^ > (2.3) 
for all u G C^(f2), s G M and all /i > small enough. 

To add the perturbation h'^"^q to the estimate (12.31) . we assume that — 1 < s < 
and use that 

Ik'^illff" J < Ik^lU^ < lklU°° II'^^IIl^ < lklU°° II'^^IIh^+i- 

To add the perturbation 

h^'^-^e^iA ■ hD)e~'^ = h^"'^\A ■ hD + iA ■ Vip) 
to the estimate (12. 3p . assuming that — 1 < s < 0, we need the following estimates 
\\{A-V^)u\\h^ < \\A-V^\\L^\\u\\fjs+i, 
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n 

\\A ■ hDu\\Hi^^ < \\hD^{A,u)\\Hs^^ + 0{h)\\{dwA)u\\H^^^ 

n 

< 0{1) J2 UAhs^ + 0{h)\\u\\H^+^ < 0{1)\\u\\hs^u 

€ scl scl scl 

When obtaining the last inequahty, we notice that by complex interpolation it 
suffices to consider the cases s = and s = — 1. 



Let 

Thus, we obtain the following Carleman estimate for a ffist order perturbation 
of the polyharmonic operator. 

Proposition 2.2. Let A e iy^'~(fi,C"), q E L°°{n,C), and Lp be a limiting 
Carleman weight for the semiclassical Laplacian on Q. Assume that m > 2. If 
— 1 < s < 0, then for h > small enough, one has 

\\C^u\\h^ > ll^illn^+i, (2.4) 

for allue C^{n). 

The formal L^-adjoint of £^ is given by 

Notice that if is a limiting Carleman weight, then so is —ip. This implies that 
the Carleman estimate (12. 4p holds also for the formal adjoint £* . 

To construct complex geometric optics solutions we need the following solvability 
result, similar to [7]. The proof is essentially well-known, and is included here 
for the convenience of the reader. In what follows, we denote by H^^y{Q) the 
semiclassical Sobolev space of order one on Q, equipped with the norm 

= Il'"lli2(f7) + II^^^IIl2(q). 

Proposition 2.3. Let A E iy^'~(fi,C"), q E L°°(fi,C), and if be a limiting 
Carleman weight for the semiclassical Laplacian on fl. Assume that m > 2. If 
h > is small enough, then for any v E L'^{Q), there is a solution u E II^{Q) of 
the equation 

C^u = V in Q, 

which satisfies 

.. .. C .. .. 
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Proof. Consider the following complex linear functional 

By the Carleman estimate (12 ■4p for the formal adjoint £*, the map L is well- 
defined. Let w G C~(fi). We have 

C 

1^(^1,^)1 = \{w,v)l2\ < \\w\\l2\\v\\l2 < —\\C*^w\\jj-i\\v\\L2, 

showing that L is bounded in the H^^-noim. Thus, by the Hahn-Banach theorem, 
we may extend L to a linear continuous functional L on H^^(R"') without increas- 
ing the norm. By the Riesz representation theorem, there exists u G H^{W^) such 
that for all w G //"^(R"), 

~ C 
L{w) = {w,u)(^H-\m), and \\u\\hi^^ < —\\v\\l2. 

Here (■, ■)(^H'^,h^) stands for the usual L^-duality. It follows that C^u = v in Q. 
This completes the proof. □ 

Our next goal is to construct complex geometric optics solutions of the equation 
'C.A,qU = in Q, i.e. solutions of the following form, 

u{x, C; h) = e~i^{a{x, () + hr{x, (; h)), 

where C e such that C ■ C = 0, |ReC| = |IniC| = 1, the amplitude a G C°°{U), 
and the remainder satisfies = ^(1)- 

Consider 

e^h^'^CA^ge-^ = {-h^A-2iC-hV)"' + h^"'-^A-hD + h'^'^-^A-C + h^"'q. (2.5) 
Since m > 2, in order to get 

e^h^"'CA,gie^a) = 0{K^+^), (2.6) 

in L^(f2), we should choose a G C°°(fi), satisfying the following first transport 
equation, 

((.V)'"a = in VL. (2.7) 

This is clearly possible. Having chosen the amplitude a in this way, we obtain 
the following equation for r, 

e^h^"'CA,ge — ^e^hr = -e^ h^"'CA,gie^ a) in n. (2.8) 

Thanks to Proposition 12.31 and (12. 6p . for /i > small enough, there exists a 
solution r G H\n) of (ESD such that ||r||j;^i^^ = 0(1). 

Summing up, we have the following result. 



8 



KRUPCHYK, LASSAS, AND UHLMANN 



Proposition 2.4. Let A e W^'°°{n, C"), q e C), and ( e be such that 

C ■ C = and I Re CI = |ImC| = 1. Then for all h > small enough, there exist 
solutions u{x, h) G H^{Q) to the equation CA,qU = in Q, of the form 

u{x, C; h) = e~h^(a(x, Q + hr{x, C; h)), 
where a(-,C) G C^iQ) satisfies (EZD and Mh^^^ = 0{1). 

Remark 2.1. In what follows, we shall need complex geometric optics solutions 
belonging to H^'^{Q). To obtain such solutions, let f2' DD Q be a bounded do- 
main with smooth boundary, and let us extend A G W^'°°{Q, C") and q G L°°{Q) 
to W^'°°{fl' jC"') and L°^(fi')-functions, respectively. By elliptic regularity, the 
complex geometric optics solutions, constructed on Q', according to Proposition 
[221 belong to H^'^iQ). 

Remark 2.2. We shall also consider the complex phases x ■ C, with ( depending 
shghtly on h, to be precise, such that ( = + ^(h) with C''°^ £ C" being 
independent of h. In this case it follows from (12.51) that we can construct complex 
geometric optics solutions with amplitudes a = a{x, C^'^^), which are independent 
of h and which satisfy the following transport equation 

(C(°^ ■ V)'"a = in Q. 

3. Proof of Theorem 11.11 

The first step is a standard reduction to a larger domain, see ^U] . 

Proposition 3.1. Let Q CC Q' be two bounded domains in R"" with smooth 
boundaries, and let A^^\ A^^^ G W^'^^Q' X"") , q^^lq'-^^ G L°^(fi',C) satisfy A^^^ = 
and = q('^ m fi' \ n. IfC%,^^,,, = C%,^^,,„ then C%,^^,,, = C%,,^^,,,. 

Proof. Let u' G H'^'^{Vl') be a solution of £^(i) = in Vl' . Since C^(i) = 
C^(2)^(2), there exists v G H'^"^{Q), solving £^(2) ^(2)^ = in and satisfying 
7f = 7m' in dQ and jv = ^u' in dQ. Setting 

f If in n, 

^ ~ [u' in n'\n, 

we get v' G H^"'{Q') and £A(2),,(2)t;' = in fi'. Thus, C%,,^^^,, C C^;2)_,(2,- The 
same argument in the other direction shows the claim. 

□ 

Recall that A^^^ G W^'°^{W, C") n£'(n, C"), j = 1, 2. Let 5(0, R) be a open ball 
in M", centered at the origin and of radius R such that Q CC -8(0, R). We extend 
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q{j) g j = 1, 2, by zero to -6(0, R) and denote these extensions by the 

same letters. According to Proposition I3.H we know that C^i)^Ji) = C^^i^J^^ . 

We shall need the following consequence of the Green's formula, see [1], 

{C-A,qU,v)L\B{0,R)) = iu,C*AqV)L^B{0,R)), U, V G H'^"' {B{0, R)) , -fU = = 0, 

(3.1) 

where C^^ = CA^i-iv A+q- 

Let Ui G i/^™(i3(0, /?)) be a solution to g(i)'Ui = in B{0,R). Then there 
exists a solution U2 G H'^"^{B{0, R)) to £^^(2) ^(2)^2 = in B{0,R) such that 
■yui = 7M2 and 7W1 = 7^2- We have 

'CA(i),,(i)(«i-«2) = (A(=^^-^^'^)-I^«2 + (g^'^-g^'V2 in B{0,R). 
Let ^; G H'^"'{B{0, R)) satisfy 

^A(i),,a)^ = in S(0,i?). 

Using (13.11) . we get 

((A^^) - A^^^) ■ Du2)vdx + [ {q^^^ -q^^^)u2vdx = Q. (3.2) 

To show the equalities A'^^'^ = A^'^'^ and g'-^^ = (/^^^ the idea is to use the identity 
(13. 2p with U2 and v being complex geometric optics solutions. To construct these 
solutions, let (^,/ii,/i2 ^ be such that = |/i2| = 1 and /^i ■ /i2 = /^i ■ ^ = 
Ai2 ■ ^ = 0. Similarly to [38], we set for > small enough. 



C2 = Y + Y 1 - /i^— /^i + «/i2, Ci = -y + Y 1 - /i^— /"I - «/^2, (3.3) 
so that (j ■ Q = 0, j = 1, 2, and C2 — Ci = 

Then by Proposition 12.41 Remarks 12.11 and 12. 2^ foTh>0 small enough, there 
exist solutions W2(x, C2; h) G H'^"'{B{Q, R)) and w(x, Ci; /i) G H'^"'{B{Q, R)), to the 
equations 

£^(2) g(2)M2 = in 5(0, i?) and £^(i)_g(i)t; = in 5(0, i?), 
respectively, of the form 

U2{xX2]h) =e~h-{a2{x, 1x1+11x2) + hr2{xX2\h)), 

^ (3-4) 
v{xX2\h)=e h (ai(x,^i -«/i2) + /iri(x, Ci;/i)), 



where the amplitudes ai(-,/ii + 2/12), a2(-, /^i — ^^^2) £ C°°(-B(0, i?)) satisfy the 
transport equations, 

((/ii + i/i2)- V)™a2(a;,//i+z//2) = 0, ((;Ui -2;U2) ■ V)"'ai(x, ^1 -2^2) = in B{0,R), 

(3.5) 
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(3.6) 



+ 



B{0,R) 
(O.R) 
B{0,R) 



(^(2) _ ^(1)) . ^e"-«(a2 + hr2)(ai + hfAdx 
h 



\\T^\Ui^ = 0{l), J = 1,2. 
Substituting U2 and f , given by fl3.4l) . in 03.21) . we get 
= 

+ f (A(^)-A«)-e"-«(Da2 + /iDr2)(ai + /ifi)rfx 
(g(2) _ g(i))e*^-«(a2 + hr2){ai + hfi)dx. 

) 

Multiplying (13. 7p by h and letting /i — )■ +0, we obtain that 

(A(2)(2;)-AW(x))e*'^■^a2(x,A^l+^/i2)al(x,/il - i/i2)rfx = 0. (3.8) 



(3.7) 



(/il+i/i2) 



B{0,R) 



Here we use (13. 6p and the fact that a-,- G C°°(S(0, i?)), j = 1, 2, to conclude that 



/ r2ridx 
Jb{q,r) 



< \\r2\\L4n\\L^ < 0{i), 



a2ridx 



B{0,R) 



<o{i), 



{hDr2)aidx 



<0(l)||/iDr2|U2<0(l). 



B{0,R) 



Substituting ai = 02 = 1 in (13. 8p . we have 



(i"i + i/U2)- / (A(') - y4('))e*"-«c/z = 0. 



(3.9) 



B(0,R) 



Similarly to [37], we conclude from (13. 9 p that 

dj{A^^^ -A^P)-dk{Af -Af'>)=0 in 5(0, i?), 1 < i, < n. (3.10) 

For the convenience of the reader, we recall the arguments of [37] . Indeed, (13. 9p 
implies that 

/i- (A(2)(0 -^(0) = forall/i,eeM", /i-^ = 0, (3.11) 

where A^^") stands for the Fourier transform of A^^\ Let ^ = (fi, . . . ,^n) and for 
j k, 1 < j,k < n, consider the vectors fi = n{^,j,k) such that fij = —C.k, 
fJ'k = other components of fj, are equal to zero. Thus, fj, ■ ^ = and 

(EUD yields that 

■ (^f (0 - - a. ■ (^(0 - ^p(O) = 0. 

This proves (I3.10p . 
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It follows from ( I3.10p that the function (/? G C^(i?(0, R)), given by 

^{x) = [ - A^^\tx)) ■ xdt, 



satisfies 

^(2)_^(i)_y^ in B{0,R). 

bmce 

^(2) _ ^(1) = in a neighborhood of the boundary dB{0,R), we conclude 
that tp is a constant, say c G C, on dB{0, R). Thus, considering ip — c, we may 
and shall assume that y9 = on dB{0, R). 

Let us now show that A^-^'^ = A^'^\ To that end, consider fl3.8p with ai(x,/ii — 
ifj,2) = 1 and 02 (x, /ii + ifi2) satisfying 

((/ii + 2/22) ■ V)a2(x,/ii + i/i2) = 1 in B{0,R). (3.12) 

Notice that such a choice is possible thanks to (13. 5p . Here we also remark that 
(I3.12P is an inhomogeneous 9-equation and we may solve it by setting 

1 f 1 

27r J^2 yi + iy2 



where x ^ C^(]R") is such that x = 1 near 5(0, R). 
We have from (13. Sp . 

(^1 + z/i2) • / (V(/j(x))e"'^a2(x,/ii + 2^2)(^a^ = 0. 

iB(0,iJ) 

Integrating by parts and using the facts that = on dB{0,R) and /ii ■ ^ = 
■ ^ = 0, we obtain that 

= / (/7(x)e*'' '^((//i + Z/X2) ■ V)a2(x, /ii + in2)dx = / ip{x)e''''^dx. 
Jb{o,r) Jb(o,r) 

Thus, = in 5(0, i?), and therefore, A^^^ = A^^). 

It is now easy to show that q^^^ = q^'^\ To this end, we substitute A^^^ = A^'^^ 
and Oi = 02 = 1 in the identity (13. 7p and obtain that 

(g(2) _ g(i))e^^-«(l + /ir2)(l + hfi)dx = 0. 

B(0,/?) 

Letting h — )■ +0, we get g'^^HO = for all ^ G M", and therefore, q^-^^ = q^"^^ 

in B{0,R). This completes the proof of Theorem II. 1[ 
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4. Boundary reconstruction and proof of Theorem 11.21 

The proof of Theorem 11.21 will be obtained by reducing it to Theorem 11.11 with 
the help of the following boundary determination result. 

Proposition 4.1. Let Q C M", n > 3, be a bounded domain with C°° -boundary , 
and let A G C°°(n,C") and q G C~(n,C). The knowledge of the set of the 
Cauchy data CA,q determines the values of A on dVt. 

When proving Proposition 14.11 it will be convenient to rewrite the equation 

CA,qU = {{-^)'^ + A-D + q)u = Q in fi, m > 2, (4.1) 



m—l 



second order system. Introducing 

ui =u, U2 = (-A)m, . . . , Um = (-A) 
we get 

{-A® I + Ri{x,D^) + Ro{x))U = in fi, 
where / is the m x m-identity matrix, U = {ui,U2, ■ ■ ■ , UmY, and 





(4.2) 



Ri{x,D^) 







... 0\ 



... 

\A{x)-D, ... 0/ 



/ -1 




Ro{x) 



A 







... -1 
\q{x) ... / 

The set of the Cauchy data for the system (14. 2 p is defined as follows, 

{{U\dn,dM\dn) ■■ U G {H^n)^, U solves (O}- 

When considering the system (14.20 near the boundary, we shall make use of the 
boundary normal coordinates. Let {x^, . . . , x^) be the boundary normal coordi- 



nates defined locally near a point at the boundary. Here x' 



ix, ...,x 



n-1 



) is 



a local coordinate for d^l and x„ > is the distance to the boundary. In these 
coordinates, the Euclidean metric has the form, see I29l. 



n-1 



.n\2 



g = gjk{x)dx^ dx^ + {dx 
j,k=i 

and the Euclidean Laplacian is given by 

-A = Dl„ + iE{x)D^n + Q(x, D^,). 

Here E G C°° and Q{x,Dxi) is a second order differential operator in D^/, de- 
pending smoothly on x"^ > 0. The principal part Q2 of Q is given by 

n-1 

j,k=i 
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Here [g^^) is the inverse of the matrix {gjk)- We remark that here E{x) and 
Q{x, D^i) are known. It follows therefore that the set of the Cauchy data for the 
system (14.21) and the set of the Cauchy data for the equation (14. ip coincide. 

In the boundary normal coordinates, the operator in the system (14. 2 p has the 
form, 

P(x, D) := {Dl„ + iE{x)Dxn + Q{x, Dx'))®I + (x, D^') + Rf^ (x)D^n + Rq{x) . 

The matrix Rq is just Rq, expressed in the boundary normal coordinates, R^i^ 
and R\ are m x m matrices all of whose entries are zero, except for the entry 
(m, 1), which is equal to Yl^Zi ^j^xi and A^, respectively. The functions Aj, 
j = 1, . . . , n, are the components of the vector field A in the boundary normal 
coordinates. 

We have the following result, which is a direct analog of [29| [33] . 

Proposition 4.2. There is a matrix-valued pseudodifferential operator B{x, Dx') 
of order one in x' depending smoothly on such that 

P{x,D) = {Dxn®I+iE{x)®I+RP{x)-iB{x,Dx')){Dxn®I+tB{x,Dx')), (4.4) 

modulo a smoothing operator. Here B{x, Dx') is unique modulo a smoothing term, 
if we require that its principal symbol is given 

Proof. Combining (14. 3 p and (14. 4p . we see that we should have 

B^{x, Dx') + i[Dxn ® /, B{x, Dx')] - E{x)B{x, Dx') + iRf\x)B{x, Dx') 

= Q{x, Dx') ®I + R']^\x, Dx') + Ro{x), 

modulo a smoothing operator. Let us write the full symbol of B{x, Dx') as follows, 

1 

j=-oo 

with bj taking values in m x m matrices with entries homogeneous of degree j in 
^' = (^1, . . . ,^„_i) e M""^ Thus, implies that 

2 ^ 11 

E ( E ^d^'b,D:,b,) - E{x) E b,{x,a+ E d.^b,{x,a (4.6) 

l=-oo j+k—\a\=l ' j=—oo j=—oo 

\a\>0, j,k<l 

1 



+iRf\x) E b,{x,0 = Q2{x,C)i + Qi{x,0i + Ri\x.0 + Mx) 

j=-oo 
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Here Qi{x, = Q{x, C,') ~Q2{x, is homogeneous of degree one in Equating 
the terms homogeneous of degree two in fl4.6p . we get 

so we should choose bi{x,^') to be the scalar matrix given by 

hix,e) = -VQ2{x,e)I- (4.7) 
Equating the terms homogeneous of degree one in (14. 6p . we have an equation for 

2hbo + dpiD:,b, - Eh + zRfh^ + d,„h = Qiix,C) + (4-8) 

|a| = l 

which has a unique solution. The terms bj, j < —1, are chosen in a similar fashion, 
by equating terms of degree of homogeneity j ' + 1 in (14. 6p . This completes the 
proof. 

□ 

Using Proposition 14.21 and the same argument as in the proof of [291 Proposition 
1.2], we conclude that the set of the Cauchy data for the system (14. 2 p determines 
the operator B{{x' ,0), D^i) modulo smoothing. It follows therefore, from fl4.7p 
and f l4.8p that the set of the Cauchy data determines the expressions, 

iA,,{x', 0) VQ2{{x', 0), e) + Yl M^'' 0)0. e'eM"-^ 

This completes the proof of Proposition 14.11 

Let r2 be a bounded domain in M" with C°°-boundary such that Q GG Q. Then 
by Proposition 14.11 we can extend A^^\ j = 1,2, to compactly supported vector 
fields in iyi^°°(fi,C") such that A^^'> = A^^) mn\n. Let us also extend q^^\ 
j = 1,2, to by zero. Theorem 11.21 follows therefore combining Proposition 13.11 
and Theorem 11.11 

5. Proof of Theorem 11.31 

Let Q C M", n > 3, be a bounded domain in R" with connected C°°-boundary. 
In this case, arguing as in Section [31 see (13.81) . we obtain that 

{l^i + ifi2)- / (A^'^^x) - A^^\x))e''''^a2{x,fii + ifi2)ai{x,fii-ifi2)dx = 0, (5.1) 
Jn 

for any fii,^2,C, £ K"' such that /ii ■ ^ = ■ = /^i ■ = and |^i| = |/i2| = 1- 
Here ai(-, /ii + 2/12), a2{-, yUi — ^/^2) G C°°{Q, C) satisfy the transport equations, 

((/ii + i//2) ■ V)™a2(a:, fii + 1^2) = and ((//i - i^2) ■ V)"'ai(x, jji - 1^2) = in fi. 
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We can go back to Section [3] and repeat the construction of complex geometric 
optics solutions, with vectors (i and C25 given by fl3.3p . where //2 is replaced by 
—^2- In this way, instead of (15. ip we arrive at 

(//i-i/X2)- / (A'^'^^x) - A^^\x))e''''^a2{x,iJi -ifi2)ai{x,idi + ii^2)dx = 0, (5.2) 
Jn 

where ai(-, — 2/12), a2{-, /^i + ^/^2) ^ C°°{Q, C) satisfy the transport equations, 

((/ii - 2/12) • V)™a2(x, /ii - i/i2) = and {{jji + i/i2) • V)"'ai{x, fii + Zy[X2) = in fi. 

Moreover, choosing 02 = ai = 1 in (15.11) and using Proposition lA.ll in the Ap- 
pendix, we get 

^(2) _ ^(1) ^ 

where (p G C^{Q). 

Choosing next ai(-,/ii — z/i2) = ai(-,/ii + ^/i2) = 1 and letting a2(-,/ii + i/U2), 
0^2 (■, /Ji — ^/U2) G C°°{Q, C) be such that 

((/ii + i^2) ■ V)a2(a;, fii + i/i2) = and ((//i - 2/12) ■ V)a2(x, /^i - z/i2) = in 
we conclude from (15.11) and (15.21) that 

(/ii ± z/i2) ■ / (Vv?)e"'^a2(x, /xi ± ifi2)dx = 0. (5.3) 

Proposition 5.1. The function ip is constant along the connected set dQ. 

Proof. Completing the orthonormal family /ii, ^2 to an orthonormal basis in M", 
/ii, /i2, /^3, • • • , ^J'n, we have for any vector x G M", 

n 

We introduce new linear coordinates in M", given by the orthogonal transforma- 
tion T : M" — > M", T{x) = y, yj = x ■ fij, j = 1, . . . ,n, and denote 

2; = Z/i+^?/2, = ^{dy,+idy^), d, = ^{dy,-idy,). 

Thus, 

(/ii z/i2) ■ V = 2^5, (/ii - 2^2) ■ V = 2d^. 
Hence, changing variables in (15. 3p . we get 

/ id-My))e'''-^giiy)dy = 0, [ idMy))e^y<g2{y)dy = 0, (5.4) 
JT{n) JT{n) 

for all ^ = (0,0, ^")> C e M'^-^ and all gi,g2 G (:7°°(Tp)) such that d^gi = 0, 
'925'2 = 0. Here ip{y) := ip{T~^y). Taking the inverse Fourier transform in (15. 4p 
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with respect to the variable we get, for all y" G ^, 

{d2^)9i{y)dyidy2 = 0, / {d^^)g2{y)dyidy2 = 0, (5.5) 

where Qy>> = T{Q) n Uy,, and Uy,, = {{yi,y2,y") : (yi,?/2) e IR^}- 

Let us show that for almost all y" G M"~^, the boundary of the set Qy" is C°°- 
smooth. To that end, consider the function / = (/i, . . . , fn-2) '■ dT{fl) — )■ M""^, 
given by {yi,y2,y") i-> y". Let p G C°°(M"',M) be a defining function of T{Q), 
i.e. ^(^]) = G M" : p(2/) < 0} and Vp{y) ^ for all y G ar(^]). Then a 
point yo G (9T(fi) is a critical point of the function / precisely when the vectors 
V/i(yo)> ■ • ■ 5 V/„_2(|/o)5 Vp(?/o) are linearly dependent. The latter holds precisely 
when {dy^p{yo),dy^p{yo)) = 0. Since fiy/ = {(?/i,|/2) : p{yi,y2,y") < 0}, it follows, 
when y" is not a critical value of /, that the boundary of Qy" is C°°-smooth. By 
Sard's theorem, the set of critical values of the function / is of measure zero and 
therefore, for almost all y" G M*^"^, the boundary of the set Qy" is C°°-smooth. 

Using Stokes' theorem, it follows from fl5.5p that for almost all y", 

/ l^gidz = 0, / ipg2dz = 0, (5.6) 



for any gi G C°°{Q,yti) holomorphic and (72 ^ C°°(f2y") antiholomorphic functions 
in Qy". In particular, taking g2 = gi in (15. 6p . we have 



(figidz = 0, and therefore, / fgidz = 0. 



Hence, 



aUy,, J any 



/ {Re(p)gdz = 0, / {lm(p)gdz = 0, 
Jany,, JdQy,, 



for any holomorphic function g G C°°{Qyii). 

Next we shall show that for almost all y", Re Lp and Im Lp are constant along the 
boundary of each connected component of Qyii. Indeed, the fact that 



{Re^)gdz = 0, (5.7) 



dfly,, 



for any holomorphic function g G C^{Qy"), implies that there exists a holomor- 
phic function F G C{^ly") such that = RetplsQ^,,. To see this, we shall 
follow the proof of [TJ Lemma 5.1], and let F be the Cauchy integral of Re ^|gn^„, 
i.e., 

F{z) = ^[ ^^^dC, zeC\dQy.. 
27^« Jan,, C- z 
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The function F is holomorphic inside and outside dVLy^ 
Sokhotski-Privalov formula states that 



hm F{z) — 



hm F{z) 



As If e C^, the Plemelj- 
Retpi^zo), zoEdQy". (5.8) 



The function C ^ (C "~ -2)^^ is holomorphic on Qyn when z ^ Qyn. Thus, ( ]5.7p 
implies that F{z) = when z Qy//. Hence, the second limit in fl5.8p is zero and 
therefore, -F is a holomorphic function on Qy// whose restriction to the boundary 
agrees with Re ^. 



Moreover, AlmF = in fi,// and ImFl 



0. Hence, F is real- valued and 



therefore, is constant on each connected component of Q,yii. Hence, Re^ is con- 
stant along the boundary of each connected component of Qyn, for almost all y". 
In the same way, one shows that Im ip is constant along the boundary of each 
connected component of Qy//, for almost all y". 

Going back to the x-coordinates, we conclude that the function if{x) is constant 
along the boundary of each connected component of the section T~^{Qyii) = 
Qr\T~^ {Uyii) , for almost all y" G M""^, where the two-dimensional plane T~^(Hy/) 
is given by 



T-\Uy») = <x = yifil + y2^2 + VjfJ'j : 2/2 



(2/3 



.,yn)- fixed 

(5.9) 

Here //i,/U2 € 1^" are arbitrary vectors such that /ii ■ /i2 = and = \fj,2\ = 1- 

Let us now prove that ip is constant along the boundary dQ. The fact that (p is 
constant along the boundary of each connected component of T~^{Qy'/), for some 
y", means that Xip = when X is a tangential vector field to dT~^{Qyii). Let 
a G dQ. Then to show that is a constant in a neighborhood of a, one has to 
establish that 

Xj(p = 0, j = l,...,n-l, 

where Xi, . . . , Xn-i is a basis of tangential vector fields to dfl in a neighborhood 
of a. 

Without loss of generality, we may assume that locally near a = (a', a^) G dfl, 



ai. 



[Xi, 



), the boundary dfl has the form x„ = 
is C°° near a'. Then for x' in a neighborhood of a', the vectors 

/ 1 \ / \ / 



Xi(x') 











,X2ix') 





V4,0(x')/ 



, . . . , Xn-l{x') 



Xn-i) where 








Va.„_,0(x')/ 

(5.10) 

form a basis of tangential vector fields to dfl in a neighborhood of a. 
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Let US show that for any Xj, j = 1, . . . ,?7, — 1, given in flS.lOp . one can find 
fii, /i2 and y" in fl5.9p . such that Xj is tangential to dT~^{Qy//). Choosing first 
fii = ei, /i2 = en, /is = e2,...,//„ = e„_i, and y" = {x2, ■ ■ ■ ,Xn-i) fixed near 
(a2, . . . ,an-i), we have that in a neighborhood of a, dT~^{Qyii) is smooth and 
is given by x„ = (j){xi,X2, ■ ■ ■ ,Xn-i)- The tangential vector field to dT~^{Qy") 
is precisely Xi{x'), and therefore, Xi(f = near a. Proceeding similarly and 
choosing /zi = e^, j = 2, . . . , n — 1, and /i2 = e„, we obtain that the vector fields 
Xj, given in flS.lOp . are tangential to the boundaries of the corresponding sections 
of Q. It follows that ip is locally constant along the boundary dQ, and since dfl 
is connected, we conclude that ip is constant along dQ. The proof is complete. 

□ 

To finish the proof of Theorem ll.3[ we argue in the same way as at the end of 
Section [31 The proof of Theorem 11.31 is complete. 

Appendix A. Characterizing curl-free vector fields 

Let C M", n > 2, be a bounded domain with boundary and let us denote 
by C^'^{Q) the space of Lipschitz continuous functions in Q and by C^'^{Q) = 
W e C\n) : Vip G CO'lp}. 

The following result, used in Section 5 in the reconstruction of the vector field 
part of the perturbation, may be of some independent interest. Notice that it 
holds without any topological assumptions on the domain Q C M". 

Proposition A.l. Let Q C M", n > 2, be a bounded domain with boundary, 
A e ^yl'~(^],C"), and 

/i ■ [ e'''-^A{x)dx = 0, for all ^, /i G M", ^ ■ fi = 0. (A.l) 
Jn 

Then there exists Lp G C^'^{Q) such that A = Vy?- 

Proof. We start by following an argument from Let xn be the characteristic 
function of Q. Then flA.ip can be written as follows, 

// • = 0, for all /X G M", ^ = 0. (A.2) 

For any vector ^ G M", we have the following decomposition, 

where RexnA^{^),lmxnA^{^) are multiples of ^, and RexnAj_{^),lmxnAj_{C) 
are orthogonal to ^. It follows from (\A.2\i that XnAj_{^) = 0, and therefore. 
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Moreover, (jAi^mplies that XnA{0) = 0. Since xnA G L^(M") n £'{n), we 
conclude that XnMO ^ n C°°(R"). Hence, a e L°°(M"). Let 

be the inverse Fourier transform of a. Then flA.SP imphes that Xn^j = DjLp in 
M"', in the sense of distribution theory, and therefore, 

Aj = Djip in Q. 

Setting tp = i-^l^, we should check that G C^'^iH). As e W^'°°{n) = 
C^'^{Q), j = 1, . . . , n, by [T9t Theorem 4.5.12 and Theorem 3.1.7] it follows that 
(f G C^'^(fi). The proof is complete. 

□ 
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